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Abstract. In this contribution to the special issue on multifcrroics we focus on 
multiferroicity driven by different forms of charge ordering. We will present the generic 
mechanisms by which charge ordering can induce ferroelectricity in magnetic systems. 
There is a number of specific classes of materials for which this is relevant. We 
will discuss in some detail (i) perovskite manganites of the type (PrCa)Mn0 3 , (ii) 
the complex and interesting situation in magnetite Fe 3 04, (Hi) strongly ferroelectric 
frustrated LuFc 2 04, (iv) an example of a quasi one-dimensional organic system. 
All these are "type-I" multifcrroics, in which ferroelectricity and magnetism have 
different origin and occur at different temperatures. In the second part of this 
article we discuss "typc-II" multifcrroics, in which ferroelectricity is completely due 
to magnetism, but with charge ordering playing important role, such as (v) the 
newly-discovered multiferroic CaaCoMnOg, (vi) possible ferroelectricity in rare earth 
perovskite nickelatcs of the type RMO3, (vii) multiferroic properties of manganites of 
the type RMn 2 05, (viii) of perovskite manganites with magnetic E-type ordering and 
(ix) of bilayer manganites. 
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1. Introduction 

Recently there has been a new surge of interest in multiferroics (MF's), single phase 
compounds in which magnetism and ferroelectrity coexist [H [2]. Such materials are 
relatively rare, which raises the fundamental question what possible mechanisms for 
multiferroic behavior can exist [3j Hj. Here we consider multiferroicity that is driven by 
different forms of charge ordering. 

1.1. Type- 1 and Type- II multiferroics 

Crudely one can divide multiferroics into two groups [5]. In type-I multiferroics 
ferroelectricity (FE) and magnetism have different origin and are often due to different 
active "subsystems" of a material. In such type-I multiferroics the magnetic order 
parameter, breaking time reversal symmetry, and the ferroelectric order parameter, 
breaking spatial inversion symmetry, coexist and have a certain coupling between them. 

In the materials that belong to the class of type-I multiferroics FE can have a 
number of possible microscopic origins [3]. A cause of ferroelectricity can be (i) the 
presence of transition metal (TM) with d° configuration, just as in BaTiOs; (ii) the 
presence of bismuth or lead where the FE is predominantly due to lone pairs of Bi 3+ 
and Pb 2+ ; or (Hi) the presence of "geometric" ferroelectricity as in YMnOs, where FE 
is caused by a rotation of rigid M-0 polyhedra (in this case MnOs trigonal bi-piramids). 
In general in these type-I materials the FE ordering temperature is much higher than 
the magnetic one. 

In type-II multiferroics ferroelectricity occurs only in the magnetically ordered state: 
ferroelectricity sets in at the same temperature as certain type of magnetic ordering 
and is driven by it. Spiral magnetic ordering, for example, can give rise to type-II 
multiferroicity [HI [7] . 

A special group of multiferroic are materials in which ferroelectricity is caused by 
the charge ordering (CO); these systems are the focus of the present paper. Most 
of the MF's and potential MF's that we will discuss are of type-I: the perovskite 
manganites (PrCa)Mn0 3 $ [2D 122] , magnetite Fe 3 4 £31 E51 EEl E3 EE], quasi one- 
dimensional organics, and the frustrated charge ordered system LuFe204 [32]. The 
complex manganites RMn 2 05 [63] are probably the first example of type-II multiferroics 
in which FE is due to the simultaneous presence of sites with different charges and with 
inequivalent bonds occurring in a magnetically ordered state. Another such example is 
the recently discovered MF in CasCoMnOe- We will discuss the observation that the 
combination of CO and magnetism may also cause MF behavior in rare earth nickelates 
RMO3 and point out that related physics explains MF behavior of perovskite manganites 
with magnetic E-type ordering as well as of particular bilayer manganites with CO . 
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Figure 1. (A) Example of a neutral one-dimensional chain that exhibiting (B) 
site-centered charge ordering, (C) bond-centered charge ordering, and (D) a linear 
combination of these two that is ferroelectric. The arrows indicate the polarization, 
which is in total zero in (B) and (C), but develops a macroscopic moment, indicated by 
the red arrow in (D). The red dashed lines in (A), (B) and (C) indicate mirror planes 
of the system. 



1.2. How charge ordering can induce ferroelectricity 

The essential mechanism by which charge ordering can lead to the appearance of 
ferroelectricity is easily explained with the help of the schematic picture shown in Fig. [TJ 
In Fig. [TjA. a homogeneous crystal (a one-dimensional chain in this case) with equal (say 
zero) charge on each site is shown. Fig. [l]B shows the same chain after a charge ordering 
in which sites become inequivalent, one set of sites has charge +e and the other -e, as 
in NaCl. This process does not breaks spatial inversion symmetry, so that the resulting 
state cannot have a net dipole moment. This is made explicit in Fig. [l]B by marking 
mirror planes of the charge ordered structure. 

Another type of charge ordering occurs when a system dimerizes, see Fig. [lp. Such 
a lattice dimerization can have different origin, e.g. a Peierls distortion. In this case the 
sites remain equivalent, but the bonds are not, as the strong and weak bonds alternate. 
One can use the terminology of a site-centered charge ordering, or site-centered charge 
density wave (S-CDW) in the case of Fig. [I]B and a bond-centered CO, or bond-centered 
charge density wave (B-CDW) in the case of Fig. [lp. Also the B-CDW structure is 
centrosymmetric and thus cannot be ferroelectric. 

If one now combines both types of charge ordering in one system, the situation 
changes drastically. The situation with simultaneous site- and bond-centered CO is 
schematically shown in Fig. [Tp. Clearly inversion symmetry is broken in this case and 
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each "molecule" (short bond in Fig. [Tp) develops a net dipole moment, so that as a 
result the whole system becomes FE. Thus solids can become ferroelectric if on top of 
site centered charge ordering also a bond dimerization occurs [8]. 

1.3. Magnetic materials in which charge ordering can induce ferroelectricity 

We will see below that the simultaneous presence of inequivalent sites and bonds can 
have a number of different origins. In some materials bonds are inequivalent just because 
of the crystallographic structure, and a spontaneous CO that occurs below a certain 
ordering temperature drives the inequivalence of the sites. Or vice versa, the material 
can contain ions with different valence, which after a structural dimerization transition 
induce FE. These two effects may also occur simultaneously, in which case there is one 
common phase transition. 

The appearance of charge ordering is in itself quite ubiquitous in transition metal 
compounds. It is often observed in systems with ions that formally have a mixed valence. 
For instance half-doped manganites like Lai^Cai^MnOs or Pr 1 / 2 Cai/2Mn03 have one 
extra electron (or hole) per two Mn's, showing charge ordering of formally Mn 3+ (d 4 ) and 
Mn 4+ (d 3 ). In these manganites CO typically extends over a large part of the doping 
phase diagram. For example in Pri_a.Caa.MnO3 CO exists for 0.3 < x < 0.85. For 
0.3 < x < 0.5 the ordering is commensurate with the same periodicity as for x = 0.5, but 
for x > 0.5 it becomes incommensurate. Magnetic order sets in at temperatures below 
the CO temperature. The idea of charge ordering developing a ferroelectric component 
in such doped manganites was first put forward in Ref. [8]. 

Iron-oxides such as magnetite (FesO^, which exhibits the famous Verwey 
transition [27J EE], and LuFe204 jl5] have a charge ordering due to the presence of both 
Fe 2+ (d 6 ) and Fe 3+ (d 5 ). In rare earth (R) nickelates of the type RMO3 the valence of 
nickel is 3+, but below a certain temperature a charge disproportionation into formally 
Ni 2 +(d 8 ) and Ni 4+ (d 6 ) takes place [291 ED]. 

Traditionally one has site-centered ordering in mind when considering CO in 
such TM compounds. On the most basic level one can view this ordering as an 
alternation of TM ions with different valencies. However bond-centered CO is also 
a real possibility. For instance, in many quasi-one- dimensional compounds with partial 
electron occupation a Peierls transition occurs, which in fact can be viewed as a bond- 
centered CO, or B-CDW [31]. Generally speaking, this possibility can be also realized 
in three-dimensional TM compounds like manganites or magnetite. There are however 
no general rules that tell us whether and when this situation occurs. Each time this 
question has to be studied separately on a microscopic level. 
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Figure 2. Different types of charge ordering in doped manganites. Left: bond- 
centered charge ordering of " Zener polarons" proposed by Daoud-Aloudine et al. [16] . 
Right: site-centered CO as proposed in the 1950's [17l 118]. 



2. Type I Multiferroics with Ferroelectricity due to Charge Ordering 

2.1. Perovskite Manganites of the type Pri^ x Ca x MnO% 

The question whether one or the other (i.e. site- or bond-centered) charge ordering 
structure (see Fig. [2j is realized in (PrCa)MnC>3 and whether a ferroelectric state can 
result from it is experimentally rather controversial. We will review this discussion after 
presenting the theoretical arguments and calculations that point into the direction of a 
possible multiferroic groundstate in these manganites. 

2.1.1. Theoretical situation In Ref. [HI E] the question of site-CO versus bond-CO was 
treated theoretically within a double exchange model that takes into account the double 
degeneracy of active e 9 -electrons, as well as the underlying magnetic structure. This 
approach had proved to be quite successful for overdoped and half-doped manganites 
before [TU1 [TT] [T2] , and could later even be extended to undoped LaMn0 3 [TS], which are 
usually considered as Mott insulators. The results of the computations in the framework 
of the degenerate double exchange model are presented in the phase diagram of Fig. [3j 
The most important result here is the blue "triangle" to the left of x ~ 0.5: in the 
simplest approximation in this part of the phase diagram the bond-centered CO is 
lower in energy than the usual site-centered one. The observation that these different 
charge ordered states are in general very close in energy is confirmed independently by 
density functional calculations [Ti] . 

A more detailed analysis shows that the situation here is actually more complicated 
and much more interesting: in this region there occurs a coexistence of site- and bond- 
centered CO. According to our general picture, such a state is FE, see right panel of 
Fig. [3j In this "triangular" part of the phase-diagram the character of the CO actually 
changes in a regular way: it evolves from a pure site-centered groundstate with magnetic 
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Figure 3. Left: Phase diagram of the degenerate double exchange model, by [5]. The 
arrow indicates a value of the ratio of the effective superexchange J and hopping t for 
which the sequence of phases as a function of doping is the experimental one. The 
letters G, 120, A, C, CE, FM and FE denote, respectively, the magnetic G, 120 degree 
Jaffet-Kittel, magnetic A, C and CE phase, the ferromagnetic (FM) and ferroelectric 
(FE) phase. Right: schematic representation of the site- and bond-centered charge 
ordered states that result in an intermediate ferroelectric state if both types of charge 
ordering are present simultaneously. 

CE-type order at x = 0.5, to an admixture of bond-centered state increasing with 
decreasing x, and finally to a pure bond-centered one at the left edge of the "triangle" 
in Fig. [3j One more point is worth mentioning. The theoretical treatment was done for 
groundstates with long-range magnetic ordering. Experimentally CO in Pro.eCacuMnOa 
sets in at ~235 K, whereas the Neel temperature is significantly lower, ~ 160 K. 
However, a recent study indicates that magnetic correlations survive in this system much 
above the temperature of a long-range three-dimensional magnetic ordering - actually 
up to T C o P3] • Such short-range magnetic correlations are sufficient for our theoretical 
approach to be valid. 

2.1.2. Experimental situation For manganites the possibility of a bond-centered charge 
ordering was first proposed in Ref. [16] . under the name of a "Zener polaron" state. 
On the basis of a detailed structural study of single crystalline Prx^Ca^MnOs (with 
x ~ 0.4) it was proposed that, contrary to the accepted charge ordering picture [T71 [18] . 
shown in the center of the right panel of Fig. |2j the low-temperature phase of this 
system is better described not site-centered, but cLS cL bond-centered CO. In this 
bond-centered CO state one extra electron is shared by the pair of neighboring Mn ions, 
moving back and forth between them and thus, by the double exchange mechanism, 
orienting the localized spins of their t2 9 electrons parallel. Hence the terminology "Zener 
polaron" - a two-site polaron with ferromagnetic coupling. The proposed bond-centered 
CO structure of Ref. p3] is schematically illustrated in the left panel of Fig. [5J 

The question whether "Zener polarons" are present in these doped manganites 
is hotly debated in the literature. Whereas on the basis of a detailed single-crystal 
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structural study the authors of Ref. [IE] claim that the bond-centered structure is 
realized at least for x = 0.4, resonant elastic X-ray scattering experiments shown that for 
x ~ 0.4 — 0.5 Mn ions are inequivalent, which was interpreted by Grenier and coworkers 
as a disproof of the Zener polaron picture [19]. A recent high-resolution transmission 
electron microscopy and electron-diffraction study, on the other hand, confirms the 
presence of a "Zener polaron" state in Pri_ x Ca x Mn03 [20] . 

From the theoretical considerations presented above it seems that this controversy 
is to an extent artificial, as it is based on an oversimplified treatment of the bond- 
centered CE structure. First of all, theoretical results show that at x — 0.5 indeed 
the conventional site-centered structure with inequivalent Mn ions is realized. But even 
for x = 0.4, where, according to our results [8j, Mn-Mn bonds become inequivalent, 
simultaneously there exists also a site-centered CO, i.e. also the Mn ions should be 
inequivalent. Thus, the theoretical results reconcile the pictures of Refs. P2] and [H] 
and show that the actual situation is in between, as it combines features of both bond- 
centered (Zener polaron) and site-centered (checkerboard type) CO. The fact that the 
actual symmetry of Pro.6Cao.4Mn03 is non-centrocymmetric Pllm [16], agrees with 
this picture and, most importantly, with the conclusions that this system should be 
multiferroic. 

A direct observation of ferroelectricity in Prx^Ca^MnOs could settle the case. 
Unfortunately the build-up of a macroscopic polarization is hindered by the finite 
conductivity of this system, also precluding direct measurement of the polarization. 
There are at the moment two experimental studies in which an anomaly (a peak) in the 
dielectric constant e was observed at the CO transition temperature [211 E2] , strongly 
suggesting the presence of ferroelectricity. Two other recent experiments were also 
interpreted as confirming the presence of FE in Pri_ x Ca x Mn03 [231 124"] . 

The anomaly of the dielectric constant at the CO transition was also found 
in a related system Pri„ z Na x Mn03 for x=0.21, which correspond to a hole doping 
n/ l =0.42 [25]. The authors argue that it is a bulk effect, and interpret it as a possible 
indication of ferroelectricity in this system. 

2.2. Fe^O^: Multiferroic Magnetite 

The manganites that we discussed in the previous section are not the only, and possibly 
even not the best example of ferroelectric and multiferroic behavior that is caused by 
charge ordering. Another intriguing and famous case is magnetite, Fe304 - the first 
magnetic material known to the mankind, see e.g. [25], with a ferrimagnetic ordering 
occurring already below ~860 K. It is also famous because of what apparently is the 
first example of an insulator-metal transition in TM oxides - the Verwey transition at 
Ty=120 K [27] [28] . a transition that in spite of almost 70 years of dedicated research 
is still not completely understood. It is much less known, however that Fe304 is 
ferroelectric in the insulating state below the Verwey temperature Ty [3H ESI ESI E3 EH] • 
So in addition to all its acclaimed fame, magnetite is probably also one of the first 
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Figure 4. Pyrochlorc lattice made by the B-sites of a spinel structure from two 
different viewpoints. In Fe3C>4 equal numbers of formally Fe 3+ and Fe 2+ ions occupy 
these sites. 

multiferroics. 

2.2.1. Charge ordering in magnetite Fe3C>4 crystallizes in an inverse cubic spinel 
structure with two distinct iron positions. The iron B sites are inside an oxygen 
octahedron and contain 2/3 of the Fe ions, with equal numbers of Fe 3+ and Fe 2+ . These 
sites by themselves form a pyrochlore lattice, consisting of a network of corner-sharing 
tetrahedra, see Fig. |4| The iron A sites contain the other one third of the Fe ions and 
are not considered relevant for the charge ordering physics. The Verwey metal-insulator 
transition at Ty=120K is apparently related to charge ordering transition of the two 
types of charges that live on the sites of the pyrochlore lattice of the iron B sites. The 
originally proposed charge ordering pattern consists of an alternation of Fe 2+ and 
Fe 3+ ions in the xy-planes at the B-sites and was later shown to be too simple. Numerous 
much more complicated CO patterns were proposed, see e.g. Ref. [32J and references 
therein. The detailed charge ordering pattern is still unresolved. 

The difficulty in determining the charge ordering structure is related to the strong 
frustration of simple bipartite ordering on a pyrochlore lattice, first discussed by 
Anderson [33] . This frustration has the consequence that there is a macroscopic number 
of charge configurations with the same groundstate energy, if only nearest neighbor 
Coulomb interactions are taken into account. This huge degeneracy remains even in 
the presence of additional constraints on the possible charge ordering patterns, the best 
known of these is the so-called Anderson constraint, where each tetrahedron is required 
to have two Fe 2+ and two Fe 3+ ions. Such a macroscopic degeneracy is expected to 
contribute significantly to the entropy of the system. It was noticed very early on 
that this situation is very similar to the frustrated order in (water) ice, in which the 
hydrogen atoms in ice are arranged and constrained according to the so-called Bernal- 
Fowler rules (39j HO] . Also in the problem of spin-ice frustration is a dominant factor, 
and its consequences have recently attracted much attention, see for instance jH] . 

It is clear that residual interactions, beyond the nearest neighbor Coulomb 
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Figure 5. Illustration of the possible origin of ferroelectricity in magnetite. 
Emphasized are the Fe chains running in [110] direction of magnetite - a pyrochlore 
lattice made by the spinel B-sites. In the xy chains there is an alternation of Fe 2+ 
and Fe 3+ ions (open and filled circles) . Simultaneously there is an alternation of short 
and long Fe-Fe bonds; shifts of Fe ions are shown by black arrows. The resulting 
polarization is indicated by the red arrows. 

repulsion, have to be present. These are important and relevant because below Ty 
magnetite is ordered into a unique groundstate. Moreover, a recent high resolution 
neutron and X-ray diffraction study by Wright and coworkers [12] concludes that the 
Anderson rule is violated in FesO^ But our main point here is that even in the very 
careful structural investigation of Wright and coworkers a centrosymmetric monoclinic 
group is used to analyze the diffraction data, which precludes the possibility of a 
ferroelectric groundstate. Thus the real structure of Fe3C>4 is apparently even more 
complicated than the ones discussed so far. 

2.2.2. Possible origin of ferroelectricity in magnetite According to Ref. (37], the 
electric polarization P of Fe3C>4 below the Verwey transition lies in the 6-direction. 
The polarization can be switched by electric field and it leads to the formation of a 
ferroelectric domain structure that can only be explained by assuming a triclinic crystal 
structure [^3] . 

The microscopic origin of ferroelectricity in magnetite remains to be unveiled, but 
one can argue (H HI] that most probably it is related to the coexistence of site-centered 
and bond-centered CO of a type that was discussed for the manganites above. Actually 
a detailed analysis of the structural data of Ref. [12] shows that the situation in this 
system may be not so far from the one illustrated in the right panel of Fig. [3j Indeed, in 
the proposed structural and charge pattern one notices that, besides the site-centered 
CO - the alternation of the formal Fe 2+ and Fe 3+ valence states - there is also a strong 
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modulation of the Fe-Fe distances (the bond lengths). In this structure, shown in Fig. [5j 
one sees that e.g. along the Fe chains running in the [110] direction (in cubic setting), i.e. 
in the xy chains, there is an alternation of Fe 2+ and Fe 3+ ions (open and filled circles 
in Fig. [5|, but simultaneously there is an alternation of short and long Fe-Fe bonds. 
This direction corresponds to the monoclinic 6-direction, in which the polarization is 
observed. In this framework each of such [110] mixed bond- and site-centered CO chains 
gives a non-zero contribution to the electrical polarization, cf. Fig. [3} 

In the analysis of the structural data in Ref. [12] the enforced centrosymmetry is 
accommodated by a certain pattern of [110] mixed chains that changes into the opposite 
pattern in the "upper part" of the unit cell that is doubled in the c-direction. In the 
enforced centrosymmetric structure Fe3C>4 is therefore antiferroelectric. How to change 
the pattern of Ref. [12] to incorporate ferroelectricity is not entirely clear at present, 
but it seems that the main structural features point very strongly in the direction of the 
mechanism proposed in Ref. [8] - the coexistence of bond-centered and charge-centered 
charge ordering - as the mechanism for ferroelectricity in magnetite. 

2.3. Frustrated and Charge Ordered LuFe20i 

Recently multiferroicity was observed in LuFe204 [1H1 US]- Despite the similarity of 
its chemical formula to the one of a spinel, the structure of LuFe204 is very different 
from a spinel. It has a double layer structure, with a triangular iron lattice within each 
layer. The average valence of Fe in it is Fe 2,5+ , similar to the B sites in Fes04, that 
also nominally have an equal number of Fe 2+ and Fe 3+ ions, or one extra electron (or 
hole) per two sites. Accordingly, charge ordering may occur in LuFe204 and indeed does 
occur at T c =330 K. 

However, because of the frustrated nature of the triangular iron lattice in each 
FeC>2 layer, the usual bipartite checkerboard Fe 2+ /Fe 3+ charge ordering is not favorable 
and another option becomes preferred. A charge redistribution between layers takes 
place, so that in each bilayer one layer, say the lower layer, has a 2:1 ratio of Fe 2+ /Fe 3+ 
and its upper counterpart has the inverse 1:2 ratio. Because of this interlayer charge 
redistribution each triangular layer can have a perfect, unfrustrated charge ordering with 
three sublattices in each triangular layer, one sublattice being occupied e.g. by Fe 3+ , and 
two others (forming a honeycomb lattice with these Fe 3+ at the center of each hexagon) 
by Fe 2+ . The charge ordering is vice versa in the other layer, see Fig. |6j As a result 
each bilayer acquires a dipole moment, shown in Fig. [6j and the total system becomes 
ferroelectric [HJ H5J H?]. As ferrimagnetic ordering sets in at T c =250 K, LuFe 2 04 is 
multiferroic below this temperature. 

Thus the FE in LuFe204 is due to a combination of two factors: the bilayer character 
of the crystal structure, and the frustrated charge ordering leading to the formation of 
charged planes (e.g., negative charging of the lower layer and positive charging of the 
upper one). As, generally speaking, in such a situation CO may be very strong, one could 
in principle expect large spontaneous polarizations - much larger than for example in 
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Figure 6. Bilayer of the FeC>2 triangular lattices in LuFe2C>4 with a schematic view of 
charge redistribution between the layers and the interlayer charge ordering that results 
in a macroscopic electric polarization indicated by the red arrows. 

type-II multiferroics with magnetically-driven ferroelectricity [2] . And indeed, the value 
of the polarization in LuFe204 is 0.24 C/m 2 - comparable to the polarization of 0.26 
C/m 2 in BaTi0 3 |3g]. 

2.4- Quasi one- dimensional organic materials 

A situation that is similar to the one discussed in the previous section - bonds 
inequivalent just because of the underlying crystal structure, and ferroelectricity caused 
by charge ordering - exists also in other materials. Probably one of the first such 
examples is the quasi-one- dimensional organic system (TMTTF^X [19]. In this material 
there is one electron (or rather one hole) per two sites (i.e. per two TMTTF molecules), 
but the molecular chains are dimerized due to a particular ordering of counter-ions 
X (X=PF 6 , AsF 6 , SbF 6 ), which gives rise to inequivalent bonds. With decreasing 
temperature at Tco ~ 50-150 K charge ordering occurs, after which the molecules (sites) 
become inequivalent, having alternating charge. The situation then becomes exactly the 
same as the one shown in Fig. [lp, and consequently the low-temperature phase of these 
materials becomes ferroelectric. This is clear from anomalies in the dielectric constant 
e, which above T c shows Curie behavior. Below Tco i n the simplest picture electrons 
(or holes) are localized at every second molecule, and the corresponding localized spins 
either undergo a spin-Peierls transition or have long-range antiferromagnetic ordering 
when temperature decreases. This behavior is seen in the generic phase diagram of 
Fig. [7| adopted from Ref. [31]. The spin-Peierls transition would effectively remove 
magnetism from the picture, but the systems with antiferromagnetic ordering should be 
classified as multiferroic. 

3. Type-II Multiferroics with Ferroelectricity due to Charge Ordering and 
Magnetostriction 

Magnetically-driven FE in type-II multiferroics may have two microscopic origins. One 
of them, probably the most common, works in systems like TbMnOa [20], Ni 3 V20g [5Tj . 
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Figure 7. Generic phase diagram of (TMTTF) 2 X or (TMTSF) 2 X molecular crystals, 
from Ref. [3XJ ■ CO - charge-ordered phase; SP - spin-Peierls one; AFM - the phase 
with antiferromagnetic ordering; SC - superconducting phase. The charge-ordered 
phase is simultaneously ferroelectric (FE) [49 . In the blue shaded region multiferroicity 
is expected to occur. 

M11WO4 [531 ES G2], in multiferroic pyroxenes [56] and in some other systems. The 
mechanism of FE in them is usually an inverse Dzyaloshinskii-Moriya effect [HH] , which 
operates in systems with non-collinear, usually spiral magnetic structures of certain 
type [61 [2] . It requires the direct action of the relativistic spin-orbit interaction. 

A second possible mechanism works also for collinear magnetic structures and does 
not require the presence of spin-orbit coupling: it is based on magnetostriction. For 
the magnetosctriction to give multiferroic behavior one usually requires the presence of 
inequivalent magnetic ions, with different charges. These, in their turn, may be either 
just different TM ions, or the same element in different valence states. This mechanism 
for multiferroicity is illustration by Fig. |8j The situation here is almost identical with 
that shown in Fig. [TJ3. Without spin ordering, the crystal structure is centrosymmetric. 
The magnetic structure itself is also inversion symmetric, but with a different inversion 
center. Taking both charge and magnetic structures together, however, the system loses 
inversion symmetry, i.e. it can become FE. The mechanism of creation of polarization 
would be the magnetostriction (MS). The MS is definitely different for ferro and antiferro 
bonds. If energy is gained by shortening of, for example, a ferromagnetic bond, the 
system will distort as shown in Fig. |8]B, and one ends up with exactly the same situation 
as shown in Fig. [Tp, with unequal bonds and different charges at opposite ends of short 
"molecules". A simple model that exhibits this behavior is the Ising model with both 
nearest-neighbor and next nearest-neighbor interactions: 

h = Ji sis z i+1 + J2 s-Si +2 . (1) 

i i 
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Figure 8. (A) Chain with alternating charges and with the spin structure 
up /up /down/ down. (B) effect of magnetostriction, shortening (in this example) the 
ferromagnetic bonds, producing a ferroelectric polarization, indicated by the red arrow 
in Fig. [If). 

If the interaction J 2 is antiferromagnetic (> 0) and sufficiently large, J 2 > l/2|Ji|, 
the magnetic ordering will be of the up-up-down-down type, see Fig. [HJ Depending 
on the sign of J 1; the bond with parallel spins (for J\ < 0) or with antiparallel ones 
(for J\ > 0) will shorten, increasing the magnetic energy gain on the corresponding 
bond. In both cases this will lead to a multiferroic state. An almost ideal realization 
of this scenario seems to be found recently in Ca 3 CoMn06 [52"] . In this system with a 
quasi one-dimensional structure the ions Co 2+ and Mn 4+ alternate along the chain, and 
magnetic structure is of up-up-down-down type. And indeed, below T/v = 16 K, Choi 
et al. observe an electric polarization in this system [52J. 

3.1. Manganites of the type RMnq,0^ 

The manganites RMn 2 05, where R is a rare earth ion [M], are the first example of type- 
II multiferroics that were discovered. In these manganites ferroelectricity is presumably 
driven by magnetostriction. It can occur because the material contains Mn sites that 
are inequivalent by virtue of the crystallographic structure of the material, with on 
top of this a magnetically driven spontaneous distortion that makes bonds unequal 
and breaks inversion symmetry. RMn 2 05 has a rather complicated crystal structure 
containing inequivalent Mn 3+ and Mn 4+ ions in different lattice positions with different 
oxygen coordination, with a magnetic ordering that breaks the inversion symmetry of 
the undistorted, non-ferroelectric lattice (65j [66] . 

Neutron and X-ray diffraction studies show that these manganites have space 
group Pbam but it is expected that in their multiferroic state the actual symmetry 
group is Pb2im, which allows for a macroscopic electric polarization along the b 
axis [H3 EH1 EH]- The orthorhombic Pbam crystal structure of HoMn 2 05 consists of 
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Figure 9. Schematic view of the crystal structure of iSVh^Os consisting of connected 
Mn 4+ 06 octahedra and Mn 3+ Os pyramids; the figure is from the review of A. B. 
Sushkov et al in this issue. The chain of Mn 3+ -Mn 3+ -Mn 4+ along 6-direction, with 
corresponding spin ordering, is shown in the upper panel. 



connected Mn 4+ C>6 octahedra and Mn 3+ C>5 pyramids (see Fig. [9]). The octahedra share 
edges and form ribbons parallel to the c axis. Adjacent ribbons are linked by pairs 
of corner-sharing pyramids. Below 38 K e.g. in HoMn 2 Os a commensurate magnetic 
structure develops with propagation vector k = (|, 0, |), and simultaneously the system 
becomes ferroelectric [70] . 

Along the 6-direction HoMn205 exhibits a charge and spin ordering that can 
schematically be denoted as a chain of Mn| + -Mn^ + -Mn| + . In the undistorted Pbam 
structure the distances d^ (between MnJ + and MnJ + ) and d^ (between Mnijj 4 " and 
MnJ + ) are the same. Ab initio calculations of the relaxed structure [721 EH E] reveal a 
shortening of distances between parallel spins Mni + and MnJ + ions - in the ferroelectric 
Pb1\m structure dff < d^, which optimizes the double exchange energy [21 El 165] . 

3.2. Origin of ferroelectricity in RMn%0^ 

There is still some controversy about the microscopic mechanism of ferroelectricity in 
this system [2]. The most plausible explanation is the one indicated above: magnetically 
induced changes in bond lengths between Mn ions with different formal valence, adding 
up to a net ferroelectric moment. The actual picture that arises from microscopic 
electronic structure calculations supports this scenario [721 (731 CI] , but is actually even 
more interesting, see below. An alternative picture is that the ferroelectricity in RM^Os 
is due to a spiral magnetic structure, like in many other type-II multiferroics [71J, 
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Figure 10. Schematic view of the two contributions to the ferroelectric polarization 
in RMJI2O5 in the uncorrelated (U=0) and strongly correlated limit (large U). In 
the latter the electronic polarization nearly cancels the ionic polarization. The labels 
"Mn 4+ / 3+ " indicate Mn ions that have a valence of more/less than 3.5+, respectively. 



but most probably the weak spiral observed in [71] is not the source, but rather the 
consequence of ferroelectricity, similar to the case of BiFeC>3. 

The electronic structure calculations show that the magnetic ordering gives rise 
to a large electronic polarization, inducing a transfer of charge from the Mn-sites to 
the Mn-O-Mn bonds, similar in spirit to the perovskite manganites that we discussed 
in a previous section. This charge transfer drives the bond-length shortening and 
gives rise to an ionic polarization as the two Mn ions involved have formally different 
charges. The transferred charge itself, however, does not induce a large dipole moment, 
because the charges accumulated on the bonds comes from opposite directions in roughly 
the same amount J72J [73j 173] . However, this changes drastically when the on-site 
Coulomb interactions between the Mn 3d electrons are accounted for [71j . It depletes 
electrons from the Mn 4+ site, making it almost closed shell t2 9 and as a consequence 
its polarization cloud disappears. The overall result is that in these strongly correlated 
multiferroic manganites a substantial electronic polarization develops, which is almost 



as large as the ionic polarization, but opposite in direction, see Fig. 10 In the end 
a small net polarization of P=82 nC/cm 2 results, in very good agreement with the 
experimental value 



3.3. Nickelates of the type RNi0 3 

In the examples of type-II MF with the magnetostriction mechanism considered above, 
the different charges of magnetic ions were determined simply by the structure of the 
compound: in the previous section we discussed different TM ions such as Co 2+ and 
Mn 4+ in CasCoMnOe or Mn 3+ and Mn 4+ in respectively the pyramids and octahedra of 
RMn 2 05. However inequivalent ions can also appear spontaneously, as a site-centered 
CO, with bonds becoming different due to magnetostriction. Possible examples of this 
are certain rare earth nickelates RNi03 (see also [2]). 

These systems are known to have a spontaneous charge disproportionation (or CO) 
below their metal- insulator transition [29| 130] . This feature is especially prominent when 
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Figure 11. Illustration of simultaneous charge disproportionation in Ni 2+ /Ni 4+ (open 
and filled circles) and magnetic ordering in RN1O3 (red and blue vertical arrows) . The 
real disproportion will be less ionic and more of the type Ni( 3-<5 ) + /Ni( 3+<5 ' + . Small 
black arrows indicate the ionic displacements along the body diagonal of the cubes due 
to magnetostriction, and the direction of the resulting macroscopic electric polarization 
P is indicated by the red diagonal arrow. 



the rare earth is small (Y, Lu, ...). The disproportionation leads to the formation of 
a rock-salt-like charge structure: alternation of, formally, Ni 2+ and Ni 4+ in consecutive 
[111] planes the cubic perovskite lattice (for simplicity we ignore orthorhombic tilting 
in this picture, which can however be important in real materials). 

It is still controversial which magnetic structure appears in these systems at 
temperatures below T N < T C o- One possibility is that it is a sequence of the same [111] 
planes ordered in fT-U fashion [57] . Another proposed structure is more complicated: 
it has the same superstructure wave vector Q=(l/4, 0, 1/4) (in orthorhombic unit cell), 
but it has Ni 2+ and Ni 4+ spins alternating in different directions, and possibly with 
non-collinear spins [29] . The second structure does not lead to FE, at least not in 
a straight forward fashion. But the first magnetic structure should naturally give a 
polarization in the same [111] direction: the distances between ferromagnetic |f and 
antiferromagnetic || planes should not be the same due to different magnetostriction 
of ferro- and antiferromagnetic bonds, which, for the alternating charges between these 



[111] layers would immediately give ferroelectricity, as illustrated in Fig. 11 



Note that this mechanism for FE is different from the one proposed in Refs. 
for manganites with E-type antiferromagnetic ordering (which is rather similar to 
the first magnetic structure of nickelates discussed above, Fig. 11). The mechanism 
presented in Ref. [59] does not require inequivalent charges of TM ions. In addition it 
would produce a polarization not in the [111] direction, but in the plane perpendicular 
to it. 
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Figure 12. (a) In-plane arrangement of Mn and atoms in orthorombic manganites 
with E-type magnetic structure. Arrows denote the direction of spins and AFM- 
coupled zigzag spin chains are highlighted by shaded areas. The Mn-O-Mn angle for 
parallel (a p ) and antiparallel (a ap ) Mn spins, large (dj and df) and small (d* and 
d 2 s ) Mn-0 bond lengths are indicated, (b) Arrows show the directions of the ionic 
displacements for Mn (left) and O (right) in the AFM-E phase. The thick arrows at 
the bottom show the direction of the resulting displacements of Mn and O sublattices 
and the polarization P [60] . 



3.4- Magnetic E-type RMnO^ Manganites 

In the perovskite manganite family RMn0 3 the magnetic E-phase was first reported for 
R=Ho [01]. The E-type magnetic structure of this multiferroic is schematically shown 



in Fig. 12 Here we wish to point out that the established picture of multiferroicity in 
E-type manganites is actually a realization of the same scenario that we have discussed 
above, of inequivalent bonds and sites combined with striction. Before explaining this 
in detail, let us retrace our steps for a moment and discuss the E-type materials in a bit 
more general context. 

As already stated, we can separate type-II multiferroics into to subcategories: 
those in which ferroelectricity is due to a particular non-collinear (spiral) magnetic 
structure [01 E] - with as microscopic driving mechanism of the polarization an inverse 
Dzyaloshinskii effect [5H] - and those in which FE is due to magnetostriction. Examples 
of the second class are e.g. RNi03 and RM112O5. These might seem very different 
systems, but they actually have, as we pointed out, the same physics as discussed 
in the main part of our paper: inequivalent sites to start with, plus bonds becoming 
inequivalent due to magnetostriction. 

In the materials discussed so far, however, one started out with TM ions with 
different charge (different valence). The suggested mechanism of FE in the E-type 
manganites [62J, supported by ab-initio calculations [60J, seems at first glance different: 
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Figure 13. (a) Schematic view of the arrangement of Mn and atoms in presence of 
a GdFeC>3 distortion. Blue arrows denote the direction of Mn spins in the "straight" 
Mn chains of the magnetic E-type ordering, cf. Fig. 12 Superexchange drives the 



Mn-O-Mn angle towards 180 degrees when spin are antiparallel and away from it 
when the spins are parallel. Green arrows indicate the resulting oxygen displacements; 
the direction of the ferroelectric polarization is indicated by the red arrow, (b) Double 
exchange also causes a ferroelectric polarization, but in the opposite direction. It drives 
the Mn-O-Mn angle away from 180 degrees when spin are antiparallel and towards it 
when spins are parallel. 



all metal ions have the same valence and charge (Mn 3+ ). However in fact the situation is 
rather similar to the one discussed in the previous section, with positively-charged Mn 3+ 
and negatively-charged 2 ~ playing the same role as TM ions with different valence is 
systems described above. This point is schematically explained in Fig. 13 which is a 



simplified representation of the real crystal and magnetic structure of Fig. 12 



As shown in Figs. 12 and [13J the oxygens are shifted away from their original 
position at the center of TM-TM bonds, due to the usual tilting of TMC>6 octahedra 
that is present in perovskites (the GdFeC>3 distortion). Note that at this point all TM- 
O-TM bonds are still equivalent. However, in the E-type magnetic structure the spin 
orientation along the TM-O-TM chains is up-up-down-down. To optimize the exchange 
in FM and AFM bonds the lattice can distort, which in this case, in contrast to Fig. [TJ 
will be achieved not by longitudinal shifts of the TM ions, but by transverse shifts of 
the oxygens, modulating the TM-O-TM angles and consequently modifying exchange 
constants. For localized electrons, according to Goodenough-Kanamori-Anderson rules 
for superexchange, the strength of a ferromagnetic bond is increased by decreasing 
the TM-O-TM angle (thus moving the oxygens in this bond down in Fig. 13 A.). To 
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strengthen an AFM bond, however, the TM-O-TM angle should increase, approaching 
180 degrees. In the magnetic E-phase this requires a shift of respective oxygens in the 
same direction, see Fig. 13 A.. 

It is interesting to note that in the case of double exchange between TM ions 
the oxygen would shift in the opposite direction: in this situation the ferromagnetic 
exchange increases with the TM-O-TM angle stretching to 180 degrees, see Fig. [13)3. 
This picture applied to the E-type manganites gives the result indicated in the right 
panel of Fig. [T2p . 

In both cases the net effect is that all negatively-charged ions (oxygens) shift in 
the same direction, which results in a net electrical polarization perpendicular to the 
chain. One can show that after such a distortion the gain in exchange energy is linear in 
displacements u, whereas the cost in elastic energy is proportional to u 2 , i.e. this process 
will always be favorable, given the crystal and magnetic structure that we discuss here. 
Clearly this is mechanism very similar to the one described earlier in this review. It 
also relies on the presence of ions with different charge (here positive TM ions and 
negative oxygens), together with bonds which become inequivalent due to a particular 
(here E-type) magnetic structure. Such bond-bending mechanism is apparently rather 
effective; the corresponding elastic modulus that counteracts such distortion is usually 
small. Therefore it can give rise to a rather large polarization 



3.5. Bilayer manganite (LaSr) 3 Mn20j 

One more example of multiferroicity that is in its nature related to the charge ordering, 
is observed in the bilayer manganite (LaSr) 3 Mn 2 07 [H3]. In this paper the evidence for 
a reorientation of orbital and magnetic stripes below 300 K is presented. After this 
happened, the situation becomes again similar to that shown in Fig. [TJ with striction 
plus CO leading to FE. Thus, the mechanism of creation of FE here is similar to the 
one discussed above, although it remains unclear so far why stripes rotate by 90 degrees 
at this transition. 



4. Conclusions 



In conclusion, we have described in this paper a generic mechanism that generates 
ferroelectricity in systems with charge ordering, which, in particular, can lead to 
multiferroic behavior in transition metal compounds. This mechanism relies on the 
simultaneous presence of inequivalent sites and bonds. This can occur spontaneously, 
due to coexistence of cite- and bond-centered charge ordering or charge density wave - 
e.g. in the case of (PrCa)Mn03. The other possibility is that one of these features is a 
property of the system itself, where either bonds or sites are inequivalent due to the very 
chemical or crystal structure of the compound. In case of inequivalent bonds, the sites 
may become different due to charge ordering; the example of that is organic compounds 
(TMTTF) 2 X or LuFe204. Such systems may be classified as type-I multiferroics - the 
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systems in which magnetism and ferroelectricity have different mechanisms and are in 
principle independent of each other. 

If the sites are different, e.g. have different valence, a bond alternation 
producing ferroelectricity may occur for example due to magnetostriction. In this case 
ferroelectricity appears only in a magnetically-ordered state with a particular magnetic 
structure: these systems are type-II multiferroics. However, in contrast to most other 
multiferroics of this class, this mechanism does not require relativistic spin-orbit coupling 
- generally a small effect in transition metal compounds. 

One of the consequences is that in principle charge ordering-induced multiferroicity 
may lead to a rather large polarization, comparable to that of good ferroelectric like 
BaTiOs and much larger than in most of the type-II (magnetically-driven) multiferroics. 
The example of LuFe204 demonstrates this. On the other hand, how strong the 
magnetoelectric coupling in this case will be and what will actually determine it, is 
still an open question. In any case, this mechanism of multiferroic behavior is certainly 
interesting and promising. Other systems of this type, besides the ones we have discussed 
in the present paper, will most probably continue to be discovered. 
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